Casimir forces between defects in one-dimensional quantum liquids 
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We discuss the effective interactions between two localized perturbations in one-dimensional (ID) 
quantum liquids. For non-interacting fermions, the interactions exhibit Friedel oscillations, giving 
rise to a RKKY-type interaction familiar from impurity spins in metals. In the interacting case, at 
low energies, a Luttinger liquid description applies. In the case of repulsive fermions, the Friedel 
oscillations of the interacting system are replaced, at long distances, by a universal Casimir-type 
interaction which depends only on the sound velocity and decays inversely with the separation. 
The Casimir-type interaction between localized perturbations embedded in a fermionic environment 
gives rise to a long range coupling between quantum dots in ultracold Fermi gases, opening a 
novel alternative to couple qubits with neutral atoms. We also briefly discuss the case of bosonic 
quantum liquids in which the interaction between weak impurities turns out to be short ranged, 
decaying exponentially on the scale of the healing length. 
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I. INTRODUCTION 

Interactions between localized defects which are medi- 
ated by the continuum they are embedded in, play an im- 
portant role in many areas of physics. Typical examples 
are the RKKY interaction between spins in a Fermi Hquid 
or the interaction between vortices in superfluids. In the 
present work, we discuss interactions between impurities 
in ID quantum Hquids. This study is motivated by the 
recent reaHzation of strongly interacting "atomic quan- 
tum wires" with ultracold gases of both bosonic % Q 
and fermionic atoms and the proposal that single 
atoms in optical traps which are embedded in a superfluid 
reservoir allow to realize an atomic analog of a quantum 
dot with a tunable coupling to the environment. Such 
quantum dots may be used to store qubits, which, un- 
der certain conditions, can be completely decoupled from 
their environment. Arrays of these dots thus appear as 
ideal candidates for quantum information processing. It 
is therefore of considerable interest to study the induced 
interactions of such dots, mediated by the environment 
they are embedded in. Similar questions arise also for 
quantum dots in solid state realizations, e.g. in carbon 
nanotubes where the interaction is mediated by elec- 
trons in the intervening wire. 

Quite generally, for both bosons and fermions, the low 
energy properties of a gapless ID quantum liquid are 
described by the so called Luttinger liquid (LL) phe- 
nomenology the effective theory is a hydrody- 
namic energy functional characterized by the velocity of 
sound u and the Luttinger interaction parameter K. In 
particular, for fermions, K = \ corresponds to the non- 
interacting case, while iiT < 1 for repulsion. For repulsive 
bosons, in turn, one has K > 1, with if — > oo in the limit 
of weak interactions, where a Gross-Pitaevskii or Bogoli- 
ubov approximation applies. As shown by Kane and 
Fisher Q (see also for a recent discussion) , the inter- 
action of a single impurity with a LL depends crucially on 



the value of K: for K > 1, the impurity is irrelevant for 
the low energy properties. A ID Bose liquid is therefore 
effectively superfluid, although there is no true conden- 
sate |3,[lH. For K < 1 the impurity changes the ground 
state of the liquid in a non-perturbative way, effectively 
cutting it into two disconnected parts. In this case, we 
will see that the induced interaction between two impu- 
rities is essentially a Casimir-like effect. Indeed, at low 
energies, two impurities at distance r define a box with 
refiecting boundary conditions for the phonon modes of 
the quantum liquid, which leads to an attractive Casimir 
interaction energy proportional to u/r with u the sound 
velocity. The case if = 1 is marginal and corresponds 
to a non-interacting Fermi gas in ID or - equivalently -, 
a system of hard core bosons, the Tonks- Girardeau gas 
In the following we will study the interactions 
mediated by the ID quantum liquid between two impu- 
rities for the various cases, including fermions with spin. 
We focus our analysis on the case of static impurities, 
while the situation of a slow time dependence, relevant 
for atomic quantum dots, where the interactions depend 
on the internal states is only discussed qualitatively at 
the end of the paper. 



II. ONE-DIMENSIONAL FERMIONIC LIQUID 

A. Non-interacting fermions 

Before considering the generic situation of impuri- 
ties embedded in a sea of interacting particles, we first 
address the marginal case K = 1 of non-interacting 
fermions. For simplicity we start from a gas of N non- 
interacting spinless fermions in the presence of two local- 
ized impurities separated by a distance r. Considering 
cold gases in atomic quantum wires, the solution of this 
problem is not just an academic exercise. Indeed, since 
fermions in a single hyperfine state have no s-wave inter- 



2 



actions due to the Pauli principle, they reaUze an ideal 
Fermi gas at sufficiently low temperatures. We assume 
that the particles are contained in a periodic box of length 
L with average density po = N/L. The (grand canonical) 
partition function of the liquid at a given temperature T 
may be expressed in terms of a functional integral over 
the Grassman fields representing the fermions: 
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Di/iD^exp(-5FL-5i). 



The corresponding action of an ideal gas is 



'FL 



/3 

da;dr 
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where r is the imaginary time running from to 1 /T = f3, 

and fi is the chemical potential (we use units such that 
h = ks = !)• The fields are anti-periodic in imaginary 
time. For short range interactions, appropriate for cold 
atoms, the interaction of the impurities with the hquid 
can be described by an additional contribution 
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(3) 



proportional to the local density at the impurity posi- 
tions Xa- Here, the index a = 1,2 labels the impurities, 
while the coupling constants Qa describe the strength of 
collisions between the atoms in the liquid and the im- 
purities. The expression for Si is based on assuming an 
effective pseudopotential for the interaction between the 
impurity and the quantum liquid. More precisely, the in- 
teraction should be replaced by a spatial integral of the 
detailed impurity potential with the microscopic density 
operator of the liquid. In the present section deahng 
with non-interacting fermions, there is no need of a high 
energy cutoff, as one can directly work with the well be- 
haved microscopic theory. However, in order to discuss 
the low energy behavior and to make contact with the 
following sections dealing with interacting fermions us- 
ing the Luttinger liquid phenomenology, we introduce 
a high energy cutoff coc- Its value can be estimated as 
u)c ~ Min{u//o, A*}) where u is the characteristic veloc- 
ity of excitations (u = vp in an ideal Fermi gas), is 
the impurity size and ji = pp/2m is the chemical poten- 
tial (at zero temperature) with pp = mvp = t^Po, the 
Fermi momentum. As will become clear from our results 
below, the coupling constants ga are then - up to a fac- 
tor vp - identical with the dimensionless backscattering 
amplitudes /i,2(7r) for fermions at the impurities. Micro- 
scopically, they are thus determined by the solution of 
the single particle scattering problem off a single impu- 
rity. In practice, an appreciable value of the backscatter- 
ing amplitude requires the impurity size to be smaller or 
of the order of the interparticle spacing, since otherwise 
the Fourier component of the potential at 2pF is close 
to zero, and hence the dimensionless couphng constants 
7q, = ga/vp vanish. Therefore, in the following, we will 
take (jJc ~ A*- 



The Grassman fields {ip, ip) are free everywhere apart 
from the points x = xi^2 and hence can be easily in- 
tegrated out by the following standard trick: first we 
formally introduce four (5-functions into the integrand 



a=l,2 

X S[tp{Xa,T) -f]a{T)]e~ 



(4) 



where (f/n , r/„ ) are the new Grassman variables describing 
the fermions at the location of the individual impurities. 
Then we introduce a set of auxiliary fields (kq,^^) us- 
ing the identity d{f) ~ / Dk cxp(« / ^/dr) to rise the 
(5-functions into the action. Finally, we integrate out the 
fermionic fields (ip^ip), which appear only quadratically, 
by Fourier transformation: 
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where 



■^0 a=l,2 



(6) 



and 



ip(xc,-X0) 



27r -icOn + 



(7) 



where = p^ /2m. — /i, and the summation is over the 
fermionic Matsubara frequencies tOn- The trivial prefac- 
tor arises from the integration over the fermionic 
fields in the absence of impurities, giving the grand par- 
tition function of the homogeneous liquid. The fields 
{KcKo) depend only on imaginary time r, or frequency 
Wn in the Fourier representation and thus the integral 
over p can be easily calculated. Since for a sufficiently 
large separation \xi — .T2I = r 'S> pj/ , the interaction en- 
ergy is small, the relevant frequencies are small compared 
to the Fermi energy uJc ^ p- ^ Vppp. An expansion to 
leading order in a;„ <C /x then gives 

q' _ \^ \^ „ „ ipF\Xa.-X„,\Sn 



X e 



-|a;„|(l-(5„ ai)lujr 



i'Y{KaVa + KaVa), (8) 



where u/r = u/r ^ p and s„ = sign(ci;„). The charac- 
teristic frequency LOr will play an important role in our 
subsequent discussions. Physically it represents the in- 
verse fiight time for a characteristic excitation in the liq- 
uid between the locations of the two impurities, which 
naturally obeys the inequality uj,. <C ojc provided the im- 
purities are much further apart than the average distance 
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between two fermions in the liquid. It is also the quanti- 
zation energy between the two impurities. In order to ob- 
tain the impurity interaction directly from the partition 
function, we integrate out the auxiliary fields (Ra^Ka)- 
This results in an action 
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fn -/ne„ \ f m 
' fnen fn 1 \ V2 



(9) 

which only depends on the four time dependent Grass- 
man fields {f]a,r]a) which describe the Fermi field at the 
impurity positions. The coefficients /„ and e„ are defined 
by 



f = ^ = ^iSrrPFr-\uJn\/uJr. 

J 71 — 95 — c 



1 - e 



(10) 



Including the contribution (EJ due to the interaction be- 
tween the impurities and the liquid, the complete expres- 
sion for the statistical sum dH) can now be written as 



-S'-S. 



(11) 



Q=l,2 



where comes from the integration of the auxiliary 
fields: 



(12) 



The total effective action S'+Si is quadratic in the Grass- 
man fields {fja,ria) and thus the integration can be done 
exactly to yield Z — Zp^Zi^Zj^, with 



n 



f3v^ 
iL 



[(/«+*7i)(/«+*72)-(/„e„)' 



(13) 

Here, the 7q, = ga/vp are the dimensionless backscat- 
tering amplitudes, characterizing the interaction of the 
impurities and the Hquid. We can now obtain the free 
energy of the ID Fermi gas in the presence of the impu- 
rities from F = — log Z/ (3 as, follows 



F=F^--Y,\og{{l-el) 

X [(/„ + j7l)(/n + 172) - (/ne„)2]}. 



(14) 



where = — logZy-^/ (3 is the free energy of the undis- 
turbed, homogeneous liquid. The expression II14II is ill 
defined as it stands, since it contains both the energy of 
zero-point fiuctuations in the g well as the formally 
divergent self-energies of the separate impurities. The 
relevant interaction energy associated with a change of 
the separation of the two impurities is given by: 

V12 = F(7„,r)-f^(0,r)-[i^(7„,oo)-F(0,oo)] 
= F(7„,r)-i^(7„,oo). (15) 



The renormalization thus requires subtracting first the 
free energy of the liquid without the impurities (7^ = 0, 
vacuum energy) and then the free energy of the system 
when the impurities are very far apart (r — > cx), self- 
energy of the impurities) . While both the vacuum energy 
and the individual self-energies are infinite in the absence 
of a cutoff, the renormalized interaction (jl5ll is finite and 
independent of the cutoff (see also the discussion below 
in section IlII A|l . 

At low temperature T ^ we can switch from sum- 
mation to integration according to dw = 27rTdn, so that 
the effective interaction energy between the impurities 
can be expressed as: 



V12 = - 



log 



1 



71726" 



-2uj I ujj--\-2ipp-r 



1 + «(7i + 72) - 7172 



(16) 



The integral can be performed analytically to yield our 
final result for the impurity interaction at T = 0: 



1/ »T- f 71726^'^"'' 

V12 = — 3?Li2 --r— 7 ; : 

2'Kr \ 1 + i(7i + 72) - 7172 



(17) 



where Li2 is the di-logarithmic function and is 
the real part. Obviously the interaction quite generally 
falls off very slowly Hke l/r with an amplitude, which 
is a strictly periodic function. Its period n/pp — p^^ 
is equal to the average inter-particle distance. This is 
a characteristic property of degenerate fermions, essen- 
tially refiecting the well known Friedel oscillations of the 
density (see below). Trivially, the interaction vanishes, if 
one of the scattering amplitudes 71,2 is zero. 

A simple expression for the renormalized interaction 
energy V12 is obtained in two limiting cases. First, if 
the interaction of the impurities with the liquid is weak 
7a ^ 1, we can expand the di-logarithm in Eq. ifTTjl to 
obtain: 



V12 = -71727; — cos(2pFf ) ■ 
Zirr 



(18) 



In the Hmit of strong impurities 7^ 1, we find in turn 
the result 



V12 = ^5JLi2 (e'^P^^) 
2iTr 



(19) 



which is completely independent of the scattering ampli- 
tudes. In this case, the interaction energy V12 can be rep- 
resented as V12 = ^f{2p^r), where f{x) = 3fJLi2(e") is 
a periodic function bounded as follows /,nin < / < /max 
where: 



max, mm 



= Li2(±l) = 



TT 



TT 

12 



(20) 



A simple way of understanding the slow 1/r-decay and 
the oscillations with period n/pp may be obtained in 
the weak scattering limit Eq. Ijl8|l . Indeed, the density 
perturbation created by a single impurity of strength 71 
at position xi is asymptotically given by: 



pi{x) w po 



7i cos{2pp\x — Xi I 
2 2t:\x — x\ I 



(21) 
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This expression for the Friedel oscillations in a spin- 
less ID non-interacting Fermi gas is valid in the limit 
where 71 ^ 1 and \x — xi\ » p^^ [l3|. Since the 
impurities couple to the local density, the interaction 
energy (excluding self-energies) of the system of two 
weak impurities is simply obtained from Ui2(ja,r) = 
52^1(3:2) + giP2{xi) where 7^ «C 1. When renormalized 
V12 = Ui2{ja,r) — i7i2(7a,oo), this interaction energy 
coincides with Eq. Ijl8|l . Alternatively, the result may be 
derived by using the random-phase approximation (RPA) 
as shown in the Appendix El 

The analysis in this section is readily generalized to the 
case of a Fermi gas with spin. In fact for non-magnetic 
impurities, such as considered in the current article, the 
two spin modes are decoupled and therefore the energy 
due to the presence of the two impurities is simply mul- 
tiplied by a factor of two. 

It should be emphasized that the calculation above, 
can be immediately extended to the case of noninteract- 
ing fermions in two or three dimensions d = 2,3, giving 
rise to an interaction energy for weak coupling of the form 



V12 ~ /i(7r)/2(7r)-^^cos(2pj.r) , (22) 

where fa (tt) are the dimensionless backscattering ampli- 
tudes of the impurities. This result follows most simply 
by considering the density fluctuations Spi{x) induced 
by a sing le impurity at position xi. As discussed, e.g., in 
Ref. [1^ they exhibit Friedel oscillations proportional to 
the dimensionless backscattering amplitude /i(7r) at the 
Fermi energy. The resulting interaction energy is then 
simply given by V12 oc /2(7r)^pi (^2). In fact, this is a 
special case of a general result that the asymptotic 
interaction between two impurities is determined by the 
product of their backscattering amplitudes. In the pres- 
ence of short range repulsive interactions between the 
fermions, we expect that the result lf22|l remains quali- 
tatively correct in the case of two and three dimensions. 
This is based on the existence of a Fermi liquid descrip- 
tion in d — 2,3, which guarantees that the low energy 
properties are qualitatively unchanged from those of a 
Fermi gas. For example, assuming that the static den- 
sity response function at 2pF is given by the particle- 
hole bubble ^3 , the renormaHzation factor Z < 1 in the 
single-particle Green function will give rise to a Fermi 
liquid correction factor in V12. This argument, how- 
ever, neglects possible vertex corrections in the density 
response which may lead to an enhancement rather than 
a suppression of the amplitude of the Friedel oscillations. 
In fact this effectively happens in the one-dimensional 
case, where the vanishing Z -factor gives rise to Friedel 
oscillations, which decay more slowly than in the nonin- 
teracting case (see below). While we are not aware of 
a quantitative calculation of the 2pp -density response in 
Fermi-liquids, it is very likely that they will give rise only 
to finite, multiplicative corrections to Eq. lf22|l . As we 
will see below, however, the situation in one dimension 



is quite different from that in d = 2,3 in the sense, that 
even qualitatively the asymptotic form of the interaction 
is not given by the Friedel oscillation picture, even for 
very weak impurities. 

Finally, we mention a recent work dealing with neutron 
matter. Bulgac et al. consider a neutron star crust, 
which is modeled as a degenerate non-interacting neutron 
gas (i.e. an ideal 3D Fermi gas) containing various kinds 
of defects or inhomogeneities (such as nuclei or bubbles) 
immersed in it. These authors compute the interaction 
energy between two defects resulting from the quantum 
fluctuations of the Fermi sea of neutrons. They obtain 
expressions similar to Eq. ll22ll , which can be interpreted 
as RKKY-like interactions between defects. In addition, 
they discuss the influence of the shape of the defects and 
consider situations with more than two defects. 



B. Spinless Fermi Luttinger liquid 

Realistic Fermi systems consist of interacting particles. 
In three and also in two dimensions, it is possible to de- 
scribe even strong interactions by Landau's Fermi liquid 
theory. As is well known, however, this concept fails in 
one dimension. Here we consider fermions with repulsive 
short-range interaction. At low-energy such a system ex- 
hibits a gapless excitation spectrum with a linear disper- 
sion, characteristic for the universality class of Luttinger 
liquids (LL) 0,ll[l3|. 

For simplicity, we start by considering spinless 
fermions, for which the low-energy description is given 
by the following hydrodynamic action: 



5ll 



2ttK 



Ax dr 



u{djf + -{drt 
u 



(23) 



Here u is the sound velocity and K the Luttinger param- 
eter. In a translationally invariant system, they obey the 
relation uK = vp Q, with vp = pp/m = ■KpQ/m the 
Fermi velocity of the associated non interacting spinless 
Fermi gas. We consider repulsive interactions for which 
K < 1. Since the Luttinger liquid description only ap- 
plies at low energies, the flelds have to be cutoff at energy 
ijJc ~ /i, where p, is the chemical potential. The associated 
cutoff length a = u/tOc is of order l/po- Of course, for a 
quantitative calculation of the scale at which the low en- 
ergy description applies, a microscopic model is needed, 
which allows to determine nonuniversal properties. For 
single impurities in Luttinger Hquids this problem has 
only recently been discussed, see ,;2fl]. Since we are con- 
cerned with the interaction at distances much longer that 
the average interparticle separation, only the low energy 
properties are relevant, which are well described by the 
hydrodynamic action II23|I . The corresponding fleld is 
related to the density of the liquid by 

p{x) « f po + ^ J [1 + 2 cos(26l + 2ppx)\ , (24) 
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where po is the equilibrium density and only the first 
harmonics are taken into account 

The interaction between the impurities and the Lut- 
tinger liquid is taken to be of the form 



rl3 

Jo , o 



(25) 



i.e. a coupling to the local density with phenomeno- 
logical scattering amplitudes (ja- Inserting the expan- 
sion H24II into this interaction, gives rise to four differ- 
ent terms. The first term is just the constant Hartree 
self-energy of the impurities, which - of course - does not 
contribute to the renormalized interaction energy V12 . In 
addition, there are terms containing dxO due to forward 
scattering. They describe quantum corrections to the 
self energies but again are irrelevant for the interaction 
V12 between two widely separated impurities. The dom- 
inant term for this interaction is the contribution pro- 
portional to cos(20 + 2pFx), which is due to backward 
scattering. In addition there are higher order terms like 
dxO cos{29 + 2ppx), however these are less relevant in the 
renormalization group (RG) sense [23|. Taking only the 
most relevant part of the interaction, the coupling be- 
tween the impurities and the LL leads to the following 
nonlinear contribution to the action 



s^[e] 



/ dr ^ 2gaPoCOs[2V7TKQ{Xa) + 2pFXa] 
•^0 a=l,2 



(26) 

wh ere w e introduced the renormalized field Q hy = 
QVttK. The complete statistical sum of the system can 
again be represented by a functional integral. To perform 
the integration over the field Q we use the same approach 
as for the ideal Fermi gas: first we introduce the new vari- 
ables 9(xi,r) — 6i(t) and 9(x2,t) = 62(t) and then 
insert the two (5-functions into the functional integral: 

z = /oe [| De„5[e(xa)-eje-^L^-5'. (27) 

" -1 10 



We then transform the S functions into the functional 
integrals over auxiliary fields and perform the Gaussian 
integration 

Z - Zl^Z^ J DeiDe2e-^"«-^'[^°l , (28) 



where the effective action for the real fields 61,2 is 

ScS = ^(©1,-ri, 62,-1 



ei,„ 
e2,„ 



(29) 

with the summation occurring over the bosonic Mat- 
subara frequencies w„, where /„ = /3|a;„|/(l — e^) and 
e„ = e"!'^"!/'^'' with = u/r. The factor Z^^ comes 
from the integration over Q{x,t) and is independent of 



the impurities, describing the homogeneous Luttinger liq- 
uid. By contrast, the factor 



Z. 



n 



°n) 



-1/2 



(30) 



which comes from the integration over the auxiliary 
fields, describes the change in the phonon modes due to 
the constraint on the Fermi fields at the positions of the 
two impurities. Similar to the noninteracting situation, 
this contribution depends on the associated characteristic 
frequency ujr and is crucial in obtaining a finite interac- 
tion energy V12 which is independent of the cutoff. 

The remaining and now non-trivial functional integral 
over the time-dependent fields Qa is of the same form 
as the one which appears in the context of quantum 
Brownian motion in a periodic potential jl^]. Indeed 
the effective action (29) basically describes two quantum 
particles subject to ohmic dissipation of dimensionless 
strength \/K which move in a periodic potential gen- 
erated by the backscattering amplitude. As has been 
shown in [23| this problem leads to a localized ground 
state \{ 1/K > 1 with small fiuctuations in the field 
6. For a quantum liquid with sufficiently strong interac- 
tions between the fermions if <C 1 and strong impurities 
7a = Qa/vF ^ 1, the functional integral Ij28ll over the 
time-dependent fields 8^ can thus be calculated using 
the stationary phase approximation (SPA): expanding 
the functions cos{2^/TTKQa + "^PpXa) from Eq. lf26ll to 
second order in the fields around one of its minimum, we 
approximate the interaction Lagrangian in the form 



5,; = 



/3 



dr 



= 1,2 



(31) 



where Ea = ^nKgaPo. Physically this means that the in- 
teraction between each of the impurities and the liquid is 
sufficiently strong to pin the local phase 6 near the value 
minimizing the potential energy. The quantities Ea play 
the role of effective frequencies for the evolution of the 
fields 9a. The approximation of the original Lagrangian 
l|26ll by the quadratic form l(3T]l is equivalent to an adi- 
abatic approximation which describes physical processes 
occurring slower than a time scale given by E~^. As the 
typical frequency of interest is ujr, the stationary phase 
approximation is valid when ^ uJr- It is thus appli- 
cable in the case of strong impurities or - equivalently - 
long distances (see below) in a strongly repulsive liquid: 
rpo% > V-ftT^ > 1. 

Within the quadratic approximation H31II . the full ef- 
fective action in Eq. Ij28|l is quadratic in 8a and hence 
can be evaluated exactly Z — Z^^^Z^Zq, where: 



= n + + /^^s) - (/„e„)'] 



-1/2 



(32) 
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The associated free energy F = — log Z/jS is given by 

n 

X [(/„ + PE^){U + - (/„e„)2]} , (33) 

where = —logZl^^/P again describes the undis- 
turbed, homogeneous Uquid. From the free energy we 
obtain the renormaHzed interaction energy between the 
two impurities V12 in precisely the same manner as in 
Eq. jnj: 

^0 V + + E2) + E1E2 

(34) 

At sufficiently low temperatures T -^ojr, the summation 
may be replaced by an integral over the real frequency 



du) log 1 



EiE2e'^'^/'^- 



- uj{Ei + E2) + E1E2 



(35) 

In the limit of long distances (or strong impurities), i.e. 
Ea S> LOr, the integral converges aX uj Ur and hence 
<C cjEa <C -Eq. For sufficiently large separations, we 
thus obtain the simple universal interaction 

Vi2^ — I dylog(l 



27rr 



-2^/^ 



24 r 



(36) 



which decays inversely with distance. The short distance 
regime, where Ea ^ tOr, can, however, not be considered 
within the stationary phase approximation (j31|l . Indeed, 
the latter is only justified if the characteristic energy scale 
<jJ ^i^r of excitations involved in the interaction does not 
exceed the effective frequencies, i.e. if E^ ^ ujr- As a 
result, for intermediate and short distances, the interac- 
tion energy V12 does not follow a simple 1/r -behavior. 
In particular, it is impossible to describe the Hmit r — > 
without properly including a cutoff or working in a mi- 
croscopic model from the beginning. It is only within 
such a more complete calculation, that the full energy 
P{.la,'r)—F{Q, r) of the two impurity problem approaches 
the self energy of the doubled single impurity case, as ex- 
pected on physical grounds. For a single scalar field in 
ID, as described by our hydrodynamic action (23), this 
calculation has recently been done by Jaffe [2^ . 

The validity of the quadratic expansion 113 111 may be 
extended to the whole relevant range K < \ oi the Lut- 
tinger parameter with the help of the self-consistent har- 
monic approximation (SCHA) [2^ 01 • In the context 
of Friedel oscillations around a single impurity in a Lut- 
tinger Hquid, this has been used by Egger and Grabert 
|25l |. It is based on making a quadratic approximation 
113 l|l for the backscattering term, however with frequen- 
cies Ea which are determined from Feynman's variational 
principle. Using Eq. 1I3TI1 as the trial action, one has to 
minimize the free energy 



\ogz + hs-s) 



Is ' 



(37) 



where S = ScS + Si and S = ScS + Si. Taking Ea as 
variational parameters, we obtain 



Ea 



AnKgaPo 



UJc_ 
Ea 



-K 



(38) 



where oJc is the high-energy cutoff. Following (23, we 
define the crossover scale tq by Ea = u/tq when both 
impurities have approximatively the same strength 71 ~ 
72 ~ 7- The SCHA is a good approximation when K < 1 
and Ea ^ ujr, i-e. at long distances r > vq. 

In the Hmit of strong impurities 7^ 3> 1, the SCHA 
frequencies Ea ~ AirKgaPo are the same as those ob- 
tained within the stationary phase approximation and 
the crossover scale is given by ropo ^ In the 

opposite limit 7a <C 1, they are given by 



Ea « AnKgaPo I 



OJc 



and the crossover scale is 



r,po-{K^r'/^'-^\poa)-^/^'-^^ 



(39) 



(40) 



where poa ^ 1 for repulsive fermions. Note the singu- 
lar behavior of the crossover scale rg ^ 00 in the limit 
if ^ 1 of non-interacting fermions. This implies that the 
regime of validity of the SCHA is moved out to extremely 
long scales rg. Quite generally, therefore, the long dis- 
tance behavior r tq of the interaction energy is al- 
ways given by the Casimir like expression l(3fijl whenever 
K < 1. The scale, however, beyond which this simple 
result applies, strongly depends on the strength of the 
backscattering amplitude and the repulsive interaction. 

In order to study the limit of weak impurities and weak 
interactions (7^ <C 1 and K < 1 close to 1) in more detail, 
we use perturbation theory. At second order in ja, and 
when r ^ p^^, we find 

V12 = -7i72(/3oa7r)^ ( - ) - — cos{2pFr)h{K) , 
\a/ Zirr 

(41) 

where h{K) = ^^^^k)^^ and T{z) is the Gamma func- 
tion. The function h{K) diverges as if ^ 1/2 and ap- 
proaches one at K = 1. The above equation was ob- 
tained under the assumption that 1/2 < if < 1. When 
if = 1, it reproduces exactly Eq. Ijl8|l obtained for 
the non-interacting Fermi gas, provided that we choose 
the following relation between the microscopic and phe- 
nomenological impurity strengths 7q = '^aPoai^- Pertur- 
bation theory breaks down when, in order of magnitude, 
the interaction energy ijlTll reaches the strong impurity 
(or long distances) result l(3fijl : 1^12! ~ Wr- This occurs 



for rpo 



7 



ipoo) ^1'^^ ^\ in agreement with 



Eq. (ESll, because 1/2 < if < 1. Therefore, the per- 
turbative result l|lT|l is valid for intermediate distances 
/Oq ^ ^ r ^ Tq. At long distances r ^ tq it is replaced 
by the Casimir-type result (j3(ill . 
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Finally, it is worth mentioning that the presence of two 
impurities implies the possibility of a tunnehng resonance 
[3,|23|. The physics of such a resonance is, however, not 
captured by the SCHA. This issue will be discussed in 
more detail in the Appendix IbI 



C. Spin- 1/2 Fermi Luttinger liquid 

In this section we generalize the above analysis by in- 
cluding the spin degree of freedom, again using a Lut- 
tinger liquid description. 

We consider N fermions in an equal mixture of spin up 



and spin down components, i.e., 7V| 



N/2. Here 



the (non- interacting) Fermi velocity is vp = pf/tti = 
7Tpo/2m with po = N/L. Taking spin into account the 
Euclidean action of the Luttinger liquid is generalized to 



E 



1 



Ax dr 



1 



(42) 

where [6^,,] is an abbreviation for [6p,9a]. This corre- 
sponds to the usual "charge" 6p and "spin" 9a fields that 
are linear combinations of the spin up and down fields: 
9p/a — -^{9^ ± 9i). The Luttinger liquid parameters 
obey the relation KpUp — vp and in addition 



K : 



Kp < 1 



and 



1, 



(43) 



where the first equation comes from considering repulsive 
interactions and the second from the SU{2) symmetry of 
the model with spin- independent interactions In 
addition K > 1/2 when considering contact interactions 
between fermions, see e.g. [ISl- The action term that 
describes the interaction with the impurities is still given 
by Eq. Ij25|l where the density is now simply the sum of 
the densities of the two spin modes that have the same 
form as Eq. lf2ijl . As stated in the previous section we 
only keep the most relevant term, which corresponds to 
backscattering by the impurity, and hence obtain: 

Si[dfj]^ 2_/ '^9°'P0 / <^T C08[V29p{Xa) +2pFXa] 
a = 1.2 -^0 

X cos[V29a{x^)] . (44) 

Now we follow the same procedure as in previous sec- 
tions and obtain 



Z = Zl^Z, / n n e-^»«[®- 



a— 1,2 fj,—p,(T 



(45) 

where we have rescaled the fields at the impurity posi- 
tions by 9 pa = \/T^KpQ^a- The effective action ^eff is 
given by 

S.s[Qpa\ =J2J2 0M5,n , (46) 



where the n refers to the Matsubara frequencies w„ and 



I^s are the elements of the matrix 



1 

— e 



--pn 



pn 



(47) 



with fpn = P\u)n\/{1 - e2„), e^„ = e-l'^"l/'^^ and Up = 
Up/r. 

In order to calculate the partition function, we again 
use the stationary phase approximation, which corre- 
sponds to expanding Si around its minima to second or- 
der in the fields, assuming that the impurities are strong, 
i.e. 7a = Qa/vp ^ 1. The nonlinear action Si is thus 
replaced by a quadratic approximation 



(48) 



where Ep^ = 2TrKpgaPQ. 

Since the charge and spin fields are now completely 
decoupled, we have Z = ZpZ^, with: 



Zp = ll[{fpn+l3Ep,){fpn + l3E^ 



p2) 



-1/2 



(49) 



The total free energy F — Fp + is then simply the 
sum of the charge and spin contribution. 

After renormalization the free energy is given by: 



/3 



fj, n>0 



1 - 



EpiEp2 e 



{ujl + Ep,){ul + Ep2) 



(50) 

At low temperature T ujp we can again replace the 
sum over Matsubara frequencies by an integral. In the 
limit of strong impurities iVp <C Epa we obtain 



(51) 



This is the straightforward generalization for spin 1/2 
fermions of the result obtained in the previous section, 
see Eq. . As discussed there, the SPA is valid only if 
LOp ^ Epa- Hence we are unable to calculate the inter- 
action at shorter distances, where ujp >• Epa- 

As in the case of spinless fermions, we can go beyond 
the SPA regime, using the SCHA. In particular we use 
S = ScS + Si and S = SeS + Si in Eq. Ij37|l : and then 
we minimize i^var with respect to Epa- In the case of 
identical impurities, i.e. 71 = 72 = 7, the values of Epa 
that minimize Fvar are such that Epi — Ep2 ^ Ep ~ 
nKpE. For large distances, i.e., for ujp ^ Ep, the SCHA 
is valid and E is given by 



E = 2gpo 1 



ttKE 



1 



ttKE 



-1/2 



(52) 



where lUc is the high energy cutoff. As in the spinless case, 
we can define the crossover scale as tq = max{up/Ep) = 
vf/ttK'^E, since K < 1. In the limit of very strong 
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impurity backscattering i?^ ^ uJc, we recover the SPA 
result, i.e. i?^ = 2TrKfj_gpo, and the crossover scale is 
ropo ^ l/K'^^. For intermediate impurity strength cj^ <C 
^ ujci we obtain 

i?, = iJ^/i^ = 2^5Poif^ , (53) 

and the crossover scale in this case is given by 

para^T^K^^iapoT^ . (54) 

where a = vpjKujc is the short distance cutoff and 
apo ^ 1- Note that, as in the spinless case, the crossover 
scale diverges rg — > oo when K ^ \. For weak interac- 
tions, therefore, the result Ij51|l is only valid at very large 
distances. 

The regime of weak impurity strength can be studied 
using perturbation theory. At second order in 7q , we find 

Fi2 =-7i72(poa^)'f-)' "^—008(2^^^^), (55) 
\a) TTr 

where h{K) ^ J^,^, (f , f ; i^; i _ ^2) 

and 2F1 is the hypergeometric function j2^. The func- 
tion h{K) is a smooth, monotonically decreasing function 
of K which diverges like h{K) ~ 2/t:K as if ^ and 
approaches one &i K = \. When K — 1, the preceding 
result reproduces exactly Eq. (tlSll (with an extra factor 
of two due to the spin degeneracy) obtained for the non- 
interacting (spinless) Fermi gas, provided that we choose 
the following relation between the microscopic and phe- 
nomenological impurity coupling constants 7q, — ^aPoai:. 
The perturbative result Ijfifill is valid for r ^ tq, with tq 
given in Eq. lf5ijl . At the crossover scale r — and for 
K <1, IV12I ~ u)r as in Eq. ll5T|l . In conclusion, therefore, 
the interaction between two impurities follows the behav- 
ior given by Eq. (I55II only for intermediate distances and 
weak interactions. By contrast, for large distances, there 
is a crossover to the universal Casimir-type interaction 
Eq. 115 111 which depends only on the velocities Up and Uo-- 

D. Discussion 

As our main result, we have shown that for separa- 
tions much larger than the interparticle spacing, the in- 
teraction energy of two impurities in a Luttinger Hquid 
of repulsive fermions (if < 1) is a Casimir-type interac- 
tion, given by a very simple universal relation, Eq. Ij36|l 
(resp. Eq. (15 111 1 for spinless (resp. spin 1/2) fermions. In 
contrast to the result Eq. ifTflll obtained for strong impu- 
rities in a non-interacting Fermi gas it does not contain 
Friedel oscillations and is independent of both the im- 
purity strengths and the interaction parameter K. The 
physical origin of this long range force is thus quite dif- 
ferent from the K = 1 case. In the non-interacting gas, 
the long range force comes from the polarization of the 



ground state. In the strongly interacting case, the Friedel 
oscillations of the ground state density around each of the 
independent impurity still exist (2^ but they are not rel- 
evant for the impurity interaction at long distances. In- 
stead the result Eq. Ij3(i|l is best understood as being the 
Casimir interaction energy of two mirrors, i.e. impene- 
trable impurities, in a phononic bath |23|. This interpre- 
tation is supported by the direct calculation of the inter- 
action energy of two mirrors in the vacuum fiuctuations 
of a ID scalar field which represents the density modes 
of the intervening quantum Hquid, see, e.g., Ref. 0. In 
fact, a similar result has previously been found for the 
force between two infinite mass beads on a string [3H . 
The Friedel oscillations are relevant for the interaction 
between two impurities only in the non-interacting case 
or at intermediate distances in the interacting Luttinger 
liquid (if < 1), see Eq. (gll and (l55)l . 

The resulting picture is consistent with the RG calcu- 
lation of Kane and Fisher for a single impurity |8| : when 
K < 1 and 7 > 0, the backscattering amplitude is renor- 
maHzed to strong coupHng in the low energy (or long 
distance r Tq) limit. The Hquid is thus effectively cut 
into pieces, with the impurities acting like perfect mir- 
rors for the acoustic modes, resulting in a Casimir force 
between them. The scale on which the impurities fiow 
to strong coupling depends on: i) the initial strength of 
the impurities 7 and ii) the fiow velocity given by 1 — if , 
see below. When the impurities are strong and the liq- 
uid is strongly interacting, the impurities fiow quickly to 
strong coupling. The associated crossover scale tq is thus 
of the order of the interparticle distance. By contrast, 
when the impurities are weak and the liquid is almost 
non-interacting, it takes very long distances to reach the 
asymptotic regime. Qualitatively, the crossover scale ro 
for two weak impurities can already be obtained from the 
scaling theory for a single impurity. Indeed, the pertur- 
bative fiow equation of Kane and Fisher gives as the 
running impurity strength: 

los^lir/a)^-'' . (56) 

For spin 1/2 fermions, the preceding equation holds pro- 
vided that the exponent 1 — if is replaced by (1 — if)/2. 
The crossover scale tq then corresponds to the distance 
at which the running impurity strength jcs becomes of 
order one, i.e. rgpo ~ -^-1/(1^^)^ in agreement with 
Eq. II4()|I . because poa ~ 1 and 1/2 < if < 1 for repulsive 
fermions with contact interactions. For longer distances, 
the impurity reaches strong coupHng and cuts the Hquid 
into disconnected pieces. 

III. ONE-DIMENSIONAL BOSE LIQUID 

In this section, we discuss the case of two identical 
impurities gS{xa) in a one-dimensional Bose liquid. In 
particular, we consider ID bosons with short range re- 
pulsive interactions gB5{x). As was first shown by Lieb 
and Liniger the dimensionless interaction parameter 
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7b = mgB/po in this problem is inversely proportional 
to the density. The strong coupling, Tonks- Girardeau, 
limit 7_B 1 is thus reached either for strong repulsion 
or at low densities. Within a low-energy effective Lut- 
tinger liquid description, the Luttinger parameter K for 
interacting bosons is larger than 1. It is related to the 
sound velocity u hy u = irpo/mK. In the weakly in- 
teracting, Gross-Pitaevskii limit 7b ^ 1, the Luttinger 
parameter diverges like K sa t^/^/tb ^ oo. For strong 
coupling 7b » 1 in turn, one finds « 1 -I- 4/7b 1. 
The singular case of non-interacting bosons (7b — 0) is 
discussed separately in Appendix O 



A. Classical ground state energy and vacuum 
fluctuations 



Before starting the explicit calculation of the renor- 
maHzed impurity interaction energy Vi2{r) in the case 
of bosons, we discuss a limitation of our quantum hy- 
drodynamic approach when treating the Bose case. We 
will shortly see, why this limitation was not discussed 
in the context of fermions. In quantum hydrodynamics, 
the ground state energy Eq is obtained as the sum of 
two terms: the classical ground state energy Eq - and the 
quantum vacuum fluctuations Eq^^\ see e.g. Ref. [3^ . 
The classical ground state energy Eq-[pq{x)] is a func- 
tional of the density profile Pq{x). Now, the presence 
of impurities or boundaries in a quantum liquid modifies 
the density profile po {x) over distances of the order of the 
healing length ^ = This in turns modifies the 

classical ground state energy in addition to the change 
in the quantum vacuum fiuctuations, which are responsi- 
ble for the Casimir-type interaction energy. In the LL 
approach, the classical ground state energy is usually 
neglected and only the fiuctuations above the classical 
ground state are taken into account. In a homogeneous 
system, the classical ground state energy is just a con- 
stant in the Hamiltonian and can therefore be safely ig- 
nored. Provided that the healing length is much smaller 
than the system size, the effect of the boundaries disap- 
pears for bulk properties. In fermionic quantum liquids, 
the healing length is of the order of the interparticle dis- 
tance 1/po- Therefore, introducing two impurities with a 
separation r ^ 1/po leaves the classical ground state en- 
ergy unchanged. However, in the case of bosons, the heal- 
ing length is much larger than the interparticle spacing, 
diverging like ~ K/po for weak interactions 7b 0. 
The condition r ^ for neglecting the contribution of 
the classical ground state energy in the calculation of the 
renormalized interaction energy between impurities V12 
thus becomes increasingly restrictive for small 7b. As we 
will see, this leaves us only with an exponentially small 
interaction, quite in contrast to the case of fermions. 



B. Weak impurities and weakly interacting bosons 

In the case of a weakly interacting Bose gas 7b ^ 1 
(corresponding to high densities) the Bogoliubov ap- 
proach is quantitatively applicable, as was shown long 
ago by Lieb and Liniger [^. The physical reason for 
that is, that the healing length ^ « l/poV7s in this limit 
is much larger than the average interparticle spacing. 
The situation is therefore essentially equivalent to the 
weak coupling limit of a 3D Bose-Einstein condensate at 
low densities, where interactions between two impurities 
have recently been discussed by Klein and Fleischhauer 
(the interaction energy V12 is called the conditional 
energy shift and denoted by A in this work). The ex- 
plicit calculation is based on the Bogoliubov approach 
and assumes that the dimensionless impurity strength 
7 = nig/irpQ is much smaller than one. Adapting the 
calculation of Ref. [IJ to a one-dimensional gas, we ob- 
tain 

2 2 

V12 il.r)^ exp {-2r/0 (57) 

to lowest order in 7. The interaction energy between 
widely separated impurities in a weakly interacting Bose 
gas thus vanishes exponentially on the scale set by the 
heahng length ^. As discussed above, this calculation as- 
sumes that the impurities modify the density profile only 
locally, which is also the physical reason for the vanishing 
interaction at r ^. 

The preceding result can again be qualitatively under- 
stood with the help of the scaling theory of Kane and 
Fisher Q for a single weak impurity in a LL. Indeed, 
when K > 1, the effective coupling of a single impurity 
is renormalized to zero in the low-energy limit. The RG 
fiow thus starts at high-energy ujc = u/a ^ p (corre- 
sponding to a short distance cutoff a ^ £_) and ends at 
a much lower energy ujr — u/r which is set by the sep- 
aration r of the two impurities. Similar to the case of 
fermions, we define a crossover scale tq by the condition 
that the effective dimensionless impurity strength at this 
scale is of order one: 

7eff - l{a/rf~^ ^ 1 ^ ro « ^71/(^-1) < ^ . (58) 

Since 7 <C 1 for weak impurities and K — tt / 00 
due to the weak interaction condition, we find rg ~ ^. 
This confirms that there is no interaction between two 
weak impurities embedded in a weakly interacting Bose 
gas when they are further apart than a distance of the 
order of the healing length. More generally, the scaling 
theory indicates that there is no long range interaction 
between weak impurities even in a strongly interacting 
Bose gas, in which K > 1. Of course the limiting case 
K ~ loi hard-core bosons is special as is the case K — oo 
of no interaction at all. The latter case is treated in 
Appendix 10 and shows, that there is no interaction be- 
tween the impurities whatever the impurity strength. In 
strong contrast to that, the limit of a Tonks-Girardeau 
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gas of hard-core bosons is equivalent to the case of non- 
interacting fermions for properties depending only on the 
modulus of the ground state wavefunction like the density 
distribution. On the basis of the calculations in section 
III AL one thus expects a long-range interaction of the form 
itTzjl between impurities in a Tonks-Girardeau gas which 
exhibits Friedel-like oscillations. In view of the fact that 
the momentum distribution of hard-core bosons is quite 
different from that of a free Fermi gas, showing no jump 
at pf, this is a quite remarkable result jl^. Another 
singular limit, where long-range interactions appear in a 
Bose liquid is that of impenetrable impurities 7 = oo. 
For arbitrary values 00 > if > 1 of the interactions, the 
interacting Bose Hquid is then cut in three disconnected 
pieces. The impurities thus act as perfect mirrors for the 
low energy phonon excitations of the intervening Bose 
liquid, giving rise to a Casimir interaction energy pre- 
cisely as in Eq. for spinless fermions. It appears, 
however, that the limit of impenetrable impurities at ar- 
bitrary energies is nonphysical, imposing strict Dirichlet 
boundary conditions on a scalar field j23|. The Casimir 
force is thus expected to be restricted to 7 = 00, while for 
any finite 7 only short range interactions should survive. 
Describing these crossovers in detail, clearly requires a 
quantitative theory of impurity interactions in ID Bose 
liquids at arbitrary values of K and 7, an interesting 
problem for further study. 



IV. EXPERIMENTAL REALIZATION AND 
DETECTION OF THE CASIMIR-LIKE FORCE 

The recent realization of one-dimensional ultracold 
Fermi gases in a strong 2D optical lattice provides 
a novel opportunity to study Luttinger Hquid effects in a 
setup with cold gases, e.g. spin-charge separation j^^]. In 
order to study whether the Casimir interactions discussed 
here might be observed in these systems, we consider an 
atomic gas of fermions in two hyperfine states. These 
two internal states play the role of (iso)spin 1/2 states. 
In principle both the sign and the strength of the inter- 
action can be controlled using scattering resonances, e.g. 
a confinement induced resonance as shown by Olshanii 
[3^ . For simplicity, we assume that the two spin states 
are equally populated = Ni — N/2. We note that 
in order to have K < 1, one needs to consider a two- 
component atomic Fermi gas. Indeed, in a single compo- 
nent Fermi gas, Pauli's principle forbids s-wave collisions 
implying that K = 1 for spinless fermions interacting via 
a contact potential. 

Following several recent ideas 0, |3 which involve 
trapping single atoms in ultracold gases, we consider an 
atomic quantum dot (AQD) Hke configuration, which 
consists of single atoms confined in a tight trap created 
either magnetically or optically, e.g. by an additional 
optical lattice. We assume that the confining potential 
can be adjusted in such a way, that it does not affect the 
atoms of the bath. The impurity atom, which is trapped 





Ib> 





Figure 1: Schematic setup of 2 AQDs coupled to a ID atomic 
reservoir. The impurity atoms (see text) in tightly confining 
potential interact with the bath when their internal level is 
\a >. Here S is the renormalized detuning and is the Rabi 
frequency coming from a laser induced coupling, see section 

El 



in a certain internal state \a) interacts with the atoms of 
the bath through s-wave collisions. In the case where two 
such AQDs are embedded in the bath and both impurity 
atoms are in state |a), the system precisely realizes the 
situation of two localized impurities interacting via a ID 
quantum liquid. Provided the Hquid consists of spin-1/2 
repulsive fermions, we expect that for distance r much 
larger than the average interparticle spacing, the inter- 
action is of the Casimir form given in Eq. Ij51|l . In princi- 
ple, using a scattering resonance may allow to reach the 
strong impurity regime 7 1 where the crossover scale 
To is even smaHer than the inter-particle distance l/po- 

A possible way to detect the interaction energy V12 (r) , 
is to do spectroscopy of a single trapped atom as a func- 
tion of the distance r to a neighboring trapped atom. In 
addition to the mean field line shifts modifying the inter- 
nal levels of the impurity atom, the Casimir interaction 
produces a Hne shift depending on distance as 1/r. For a 
quantitative estimate of this effect, we compute the en- 
ergy Ij51ll for the experimental situation realized in Ref. 
3^ There, about N - 100 ^"K atoms (per tube) form 
an atomic wire of length L ^ 10 ^m. The temperature 
can be as low as T ~ 50 nK, which is about one tenth 
of the Fermi temperature. The Fermi velocity is of or- 
der vf ^ 2.10~^ m/s and we take + Up ^ 2vf- As 
the tube length is of the order of 10 ^m and the inter- 
particle distance 1/po ~ 0.1 ^m, we assume the inter- 
impurity distance to be r ~ 1 /im, which is larger than 
the crossover length for strong impurities. This gives a 
Casimir-related line shift of the order of 1 kHz, which is 
in an experimentally accessible range. With the param- 
eters given above, the characteristic frequency uJr is of 
the order of ~ ST, which is not much larger than T as 
required for the validity of the zero temperature limit in 
which the Casimir force is obtained. 
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V. CONCLUSION 

In conclusion, we studied the long-range interaction be- 
tween two impurities mediated by the ID quantum liquid 
in which they are embedded. We found that for repulsive 
fermions, the impurities interact via a RKKY-like inter- 
action at intermediate distances and via a Casimir-Hke 
force at large distances. The crossover scale separating 
these two regimes depends on the strength of the im- 
purities and on the interactions between fermions. We 
proposed an experimental realization of such a system 
with atomic quantum dots in an ultra-cold atomic gas 
and suggested a way to detect the Casimir-type interac- 
tion by spectroscopy of a single atom in an AQD. 

An issue which is still open is to understand the inter- 
action between impurities in a strongly interacting Bose 
liquid. In particular, how does the short-range inter- 
action (on the scale of the healing length) turn into a 
long-range interactions featuring Friedel oscillations in 
the Tonks-Girardeau limit? Another issue is to assess 
the validity of the self-consistent harmonic approxima- 
tion (SCHA) used to discuss spin 1/2 fermions in section 
III CI Indeed, it is not obvious that the variational ansatz, 
which assumes decoupling of the charge and spin modes 
in presence of the impurities, is a valid starting point. 

In this paper, we studied static impurities. The situ- 
ation becomes even more interesting if one has dynamic 
impurities, as in an AQD. First, we discuss the possibility 
of internal dynamics for the AQD. We have seen that the 
characteristic frequencies of vacuum modes (excitations) 
responsible for the long-range interactions between AQDs 
are limited by uir- This means that the effective interac- 
tion potential between static impurities can be used (in 
an adiabatic approximation) for time-dependent impu- 
rity strengths 7a (i) provided that the interaction prop- 
erties change slowly compared with the time scale w^^. 
Consider, for instance, the configuration with two AQDs 
described previously, see Figure (a similar scheme for 
two impurities in a 3D Bose-Einstein condensate is dis- 
cussed in Ref. ^j). The two- level impurity atoms can 
be described as isospins 1/2 or qubits. A laser can drive 
transitions (equivalent to single-qubit gate) between the 
two internal levels \a) and \b). In the adiabatic approx- 
imation, the AQD variables are slow and can be taken 
out of the integrals so that our previous treatment to 
calculate the interaction between two impurities applies. 
Thus one can easily write an effective Hamiltonian for 
the 2 AQDs in the form 

a = l,2 ^ ^ 

+ iFi2(aW + l)(a(2) + l), (59) 

where S is the (renormalized) detuning, Q is the (effec- 
tive) Rabi frequency coming from the laser induced cou- 
pling and ai"\ (Ty"\ af^ are the Pauli matrices de- 
scribing the isospin 1/2 of each AQD a = 1,2. The long- 



range potential V12 depends, as we have seen, strongly 
on the characteristics of the bath and on the distance 
between the AQDs. In addition, the case of impurities 
with internal dynamics embedded in a spin 1/2 fermionic 
liquid is of course also relevant for discussing the RKKY 
interaction in Luttinger liquids, as discussed perturba- 
tively in fs^. 

It is also possible to imagine external motion or dy- 
namics for the impurities. The argument of adiabaticity 
also holds in the case where the distance between the 
impurities changes sufficiently slowly. This means that 
the Casimir-like interaction could be used, for example, 
to create long range attractive forces in mixtures of ID 
fermionic gases. However, when the external dynamics of 
the impurities are taken into account on the same foot- 
ing as the bath dynamics other effects could reduce, if 
not wash out, the Casimir force 0|- 
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Appendix A: RPA FOR WEAK IMPURITIES IN 
AN IDEAL FERMI GAS 

In this appendix, we give another derivation of the 
weak coupling result ifTsjl using the random phase ap- 
proximation (RPA) . The interaction between the two im- 
purities may be written as the sum of two contributions: 
a direct interaction (which is zero in the present case 
due to the short range nature of the impurities), and an 
indirect interaction induced via the polarization of the 
medium. The polarization operator in the Fourier rep- 
resentation n(u;„, q) has a singularity at g = 2pi?, which 
leads to the appearance of the long range force. Specifi- 
cally, we calculate the interaction energy of the two im- 
purities using the RPA, which works well whenever only 
states with very few particle-hole pairs are excited by the 
perturbation, i.e. in the limit 7a ^ 1. For an ideal (spin- 
less) Fermi gas the RPA Lagrangian can be represented 
as (r = and here we use the real time formalism from 
the start): 

= i /■ dtE5^(<^;,0p, ~ c.^,4<^,,), (Al) 
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where pk = X]p(2wpfe)"'^/^0pfe is the RPA expression for 
the fermion density, ojpk = {p + /2m — /2m is the 
energy of the electron-hole pair, and the summation over 
the momentum p is Hmited by the conditions: \p\ < pf 
and |p + A:| > pp. The interaction of the AQDs with the 
liquid can be written as 



S 



RPA 



2 ^ 

pk 



Vk + ^IkV^k), (A2) 



where Vk = gi + .g2exp(ifcr). The total RPA action 
^RPA_|_^RPA quadratic, and hence the fields (j)pk can be 
integrated out, so that the effective interaction between 
the impurities is given by V12 = — ^ J2k l^fepn(a; = 0, k), 
where the polarization operator is: 



n(o,fc) 



dp 



2pk + k"^ Trk 



log 



2Pf 



k — 2pp 



(A3) 

Substituting Vk, performing the integration by trans- 
forming the integral along the real k axis into the inte- 
grals along the branch cut corresponding to the singular- 
ity of the log-function, and removing the self energies of 
the separated impurities (renormalizing the interaction) , 
we find 



V12 = -I1I2-Z — cos{2pFr) 
Zirr 



(A4) 



in agreement with Eq. II18|I . 



Appendix B: COULOMB BLOCKADE AND 
RESONANT TUNNELING 



The present paper discusses the interaction energy be- 
tween two impurities in an atomic quantum wire (i.e., 
a ID quantum liquid made of atoms). In this appendix, 
we wish to make contact with related subjects in the field 
of mesoscopic conductors (or solid-state quantum wires) : 
namely Coulomb blockade and resonant tunneling of elec- 
trons in a quantum wire with two tunnel junctions or 
barriers (see, e.g., 0, 0, H^). For simplicity, we only 
consider the case of spinless fermions. 

We first consider an atomic quantum wire with 
two impurities. When the impurities are strong 
(or very far apart), the atomic quantum wire is 
cut in three disconnected pieces and we may picture 
the low energy behavior of the system as the fol- 
lowing structure: left wire/impurity /island (or central 
wire) /impurity /right wire. The Casimir effect occurring 
in this system is directly related to the energy cost of 
transferring a supplementary particle from one of the 
wires to the island. In our case (neutral atoms inter- 
acting via a short-range potential), this energy cost is a 
finite size energy ji^l equal to: 



TTVp /2K'^r 



(Bl) 



With the help of the approximate relation K k, \ -\- 
gj/nvp this energy cost can be seen as the sum of 



two contributions: a kinetic energy cost 'iTVp/2r and an 
interaction energy cost, resulting from the local interac- 
tion with the other particles on the island, g//2r, where 
gf is the forward scattering coupling constant (in the 
standard notation of the g-ology, gj — g2 = <74, see [lfl| . 
e.g.). Note that the finite size energy is of the order of 
the zero-point kinetic energy of a phonon on the island 
TTVp /2Kr. 

In the case of electrons in a solid-state quantum wire 
with two tunnel junctions (corresponding to the struc- 
ture: left electrode/barrier/island/barrier/right elec- 
trode), the Coulomb blockade is due to the energy cost 
of transferring a single electron from an electrode to the 
island ^ ||22j. However, this energy cost is not only 
due to the finite size energy Trvp/2K'^r (with K^^ w 
1 -I- e'^/eirvp, where e is the electron charge and e is an 
appropriate dielectric constant) but also gets an addi- 
tional contribution from the charging energy e^/2C of 
the capacitors (i.e. the two tunnel junctions), where C is 
the sum of the capacitance of each tunnel junction [23 |. 
The total energy cost is equal to ttvf /2K'^r+e'^ /2C . The 
charging energy is due to the long-range part of the in- 
teraction between electrons and therefore does not arise 
in the case of cold atoms interacting via a short-range 
potential. 

Another subject of comparison between the atomic 
quantum wire and the solid-state quantum wire is the 
possibility of tunne ling resonances across the double bar- 
rier structure 0, l2fil | (see for review). Indeed, 
in an atomic quantum wire with two impurities, when 
cos{pFr) =0 0, ll^ a tunnehng resonance occur: the 
energy cost to add a particle on the island vanishes and 
particles can therefore tunnel trough the impurities. The 
liquid is no more cut into disconnected pieces and we 
do not expect a Casimir effect to occur. Tunnehng res- 
onances do not appear in our calculations because they 
are not captured by the self-consistent harmonic approx- 
imation we used, as discussed in Ref. [2^ for example. 
However such resonances are infinitely sharp at zero tem- 
perature 9, 26] and therefore they do not play a major 
role and should be easy to avoid experimentally. 



Appendix C: IMPURITIES IN AN IDEAL BOSE 
GAS 



The ideal Bose gas (75 = 0) is a singular case: the 
healing length ^ diverges and boundaries are therefore 
felt over macroscopic distances. Here using quantum hy- 
drodynamics makes no sense, but of course one can ex- 
actly solve the problem from first principles (Schrodinger 
equation). At zero temperature, all bosons are in the 
same single-particle wavefunction tpo (the corresponding 
many-body wavefunction is just a product or Hartree 
state), which is the ground state of the single-particle 
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Schrodinger equation 

-T!-dlMx)+9[S{x'r/2)+6{x + r/2)]Mx) = eoM^), 
Zm 

(CI) 

where g is the impurity coupHng constant. We assume 
that the particles are on a ring of length L. The ground 
state energy Eo for N bosons is given by A^eo- For 5 > 0, 
it is bounded as follows 

eo(g = 0,r)< eo(.9, r) < eo(<? = 00, r), (C2) 

which is a direct consequence of the Schrodinger equa- 
tion. In the following we will show that both bounds 
are going to zero in the thermodynamic limit, implying 
that Eo(g,r) = for all <? > 0. Therefore, the interac- 
tion energy ^12 (5, r) = Eo{g, r) — Eo{g, r 00) vanishes 
whatever the distance between the impurities. 

On the one hand, when the impurity strength is zero 
(7 = 0, the ground state of the single-particle Schrodinger 
equation with periodic boundary conditions is just the 
constant wavefunction, which has zero energy eo{r,g = 
0) — 0. On the other hand, when g — go, the wave- 
function has to vanish on the location of the impurities, 
implying some bending of the wavefunction and a corre- 
sponding cost in kinetic energy. The ground state wave- 



function is 

^o(^) = i^raS^n(^^-:^^ \{\x\>r/2 

= if \x\ < r/2 (C3) 

where ipm = ^o{x — ±L/2) is the maximum value of the 
wavefunction and the energy is: 

^oi9 = oo,r) = ±[^y (C4) 

This quantity vanishes in the thermodynamic limit {L 
00 at fixed density po = N/L) such that L'^ r. There- 
fore eo(g = 00, r) = in the thermodynamic limit, for all 
r such that r L. Of course, this conclusion does not 
hold for an ideal Fermi gas, because fermions have to oc- 
cupy different single-particle states (following the Pauli 
principle) and therefore the average energy per particle 
does not vanish in the thermodynamic limit. 

In conclusion, there is no interaction energy between 
two impurities in an ideal Bose gas, provided that the 
distance r is much smaller than the ring size L, which is 
always satisfied in the thermodynamic limit. 
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Such a finite size energy was first discussed by Haldane 
Q in a slightly different context. He noticed that for a 
finite size Luttinger liquid, phonons are not the only type 
of low-energy excitations, but that there are also two so- 
called "zero modes". One of the zero modes corresponds to 
adding a particle to the system without adding current. 
For a system of size r, the corresponding energy cost is 
7rwjv/2r (where vn ~ u/K — vf/K'^ is a characteristic 
velocity depending on the compressibility of the system, 
which was introduced in Ref. |3|). This is precisely the 
energy cost of transferring a particle from one of the wires 
to the island, see Eq. |IbT)| . 



